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Abstract 

Bandwidth-sharing networks as introduced by Massoulie & Roberts (1998) model the 
dynamic interaction among an evolving population of elastic flows competing for several 
links. With policies based on optimization procedures, such models are of interest both from 
a Queueing Theory and Operations Research perspective. 

In the present paper, we focus on bandwidth-sharing networks with capacities and ar- 
rival rates of a large order of magnitude compared to transfer rates of individual flows. This 
regime is standard in practice. In particular, we extend previous work by Reed & Zwart 
(2010) on fluid approximations for such networks: we allow interarrival times, flow sizes and 
patience times to be generally distributed, rather than exponentially distributed. We also de- 
velop polynomial-time computable fixed-point approximations for stationary distributions of 
bandwidth-sharing networks, and suggest new techniques for deriving these types of results. 

Keywords: bandwidth-sharing, rate constraints, impatience, large capacity scaling, fluid 
limits, fixed-point approximations. 

MSC2010: Primary 60K25, 60K30, 60F17, 60G57; Secondary 90B15, 90B22. 



1 Introduction 

Bandwidth-sharing policies as introduced by Massoulie & Roberts |25 ( l21 j dynamically distribute 
network resources among a changing population of users. Processor sharing is an example of such 
a policy and assumes a single resource. Bandwidth-sharing networks are of great research and 
practical interest. Along with the basic application in telecommunications media, e.g. Internet 
congestion control, they also have recently been suggested as a tool in analyzing problems in 
road traffic 1191. 



The main issues in bandwidth-sharing related research are stability conditions and performance 
evaluation. A variety of results regarding the first topic may be found in De Veciana et al. 
[271 El], Bonald & Massoulie [5], Mo & Walrand [22], Massoulie [20], Bramson |7J, Gromoll & 
Williams [12], and Chiang et al. [8]. As for the second topic, for special combinations of network 
topologies and bandwidth-sharing policies, the network stationary distribution may be shown 
to be of a product form insensitive to the flow size distribution, see Bonald et al. [6]. However, 
in general, approximation methods must be used, which is the subject matter of the present 
paper. Fundamental papers on fluid limit approximations for bandwidth sharing-networks are 
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Kelly & Williams |18j and Gromoll & Williams [13] . some more results on fluid and diffusion 
approximations are to be found in Borst et al. [9J 0] , Kang et al. [E] and Ye & Yao [29j EQ] • The 
latter works ignore the impact of individual peak rate limitations though, as has been pointed 
out by Roberts [23] . 

To the best of our knowledge, Ayesta & Mandjes [2J were the first to deal with fluid and diffusion 
approximations of bandwidth-sharing networks with peak rate limitations. They consider two 
specific settings first without rate constraints, and then they truncate the capacity constraints 
at the rate maxima. Reed & Zwart [23] develop a different approach in the context of general 
bandwidth-sharing networks. They incorporate the rate constraints into the network utility 
maximization procedure that defines bandwidth allocations. Thus, users operating below the 
maximal rate are allowed to take up the bandwidth that is not used by other rate constrained 
users, and bandwidth allocations are Pareto optimal. Another interesting feature of this work 
is the scaling regime. In contrast to the papers mentioned above, which mostly focus on the 
large-time properties of networks with fixed-order parameters, Reed & Zwart view networks 
on a fixed-time scale letting arrival rates and capacities grow large. This large capacity scaling 
reflects the fact that overall network capacity and individual user rate constraints may be of 
different orders of magnitude. For example, it is common that Internet providers set download 
speed limitations for individual users which are typically measured in megabits per second, while 
network capacities are measured in gigabits or terabits per second. 

The present paper builds upon [23] by relaxing its stochastic assumptions: we assume general 
distribution for interarrival times and general joint distribution for the size and patience time 
of a flow (in particular, the flow size and patience time are allowed to be dependent), while [23] 
assumes a Markovian setting with independent arrivals, flow sizes and patience times. We study 
the behavior of bandwidth-sharing networks in terms of measure- valued processes that are called 
state descriptors and that keep track of residual flow sizes and residual patience times. The first 
main result of the paper is a fluid limit theorem (it generalizes the fluid limit result of [23] to 
non-Markovian stochastic assumptions). We propose a fluid model, or a formal deterministic 
approximation of the stochastic bandwidth-sharing model, and show that the scaled state de- 
scriptors are tight with all weak limit points a.s. solving the fluid model equation. We provide a 
sufficient condition for the fluid model to have a unique solution, which converts tightness of the 
scaled state descriptors into convergence to this fluid model solution. In the sense of techniques 
used in the proofs, this part of the paper is closely related to previous work on bandwidth- 
sharing [13] , processor-sharing with impatience [11] , and bandwidth-sharing in overload [3J [9] . 
The rate constraints play a crucial role in adopting these techniques. For example, the proof 
of convergence to fluid model solutions in [11] requires an additional assumption of overload to 
eliminate problems at zero. However, in our case, due to the rate constraints, the network never 
empties, and the load conditions become irrelevant. 

Our second main result, which is a new type of result for bandwidth-sharing networks, is con- 
vergence of the scaled network stationary distribution to the fixed point of the fluid model, 
provided the fixed point is unique. There is a similar result by Kang & Ramanan [IT] for a call 
center model, but the techniques of [17] are different than ours. Applying the approach of Borst 
et al. [3], we prove that in many cases the fixed point can be found by solving an optimization 
problem with a strictly concave objection function and a polyhedral constraint set, and thus is 
unique and computable in polynomial time. We also construct an example with multiple fixed 
points, which is a feature that is distinctive from earlier cited works. Besides proving new results 
for the particular model of bandwidth-sharing, we also suggest new ideas and believe that they 
can be adjusted to other models, too. In particular, we derive equations for asymptotic bounds 
for fluid model solutions (see Theorem [3]) that can be solved for a wide class of networks, and 
then asymptotic stability of the fixed point can be shown. Another interesting idea is that, in 
the stationary regime, the properties of a network depend on newly arriving flows only, since 
all initial flows are gone after some point (see Lemma [2]). Throughout this part of the paper, 
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we assume Poisson arrivals, since that guarantees existence of a unique stationary distribution. 
Poisson arrivals also imply M/G/oo bounds that are exploited heavily in the proofs. 

The structure of the paper is as follows. Section [2] describes the stochastic bandwidth-sharing 
model, and Section[3]introduces its deterministic analogue, the fluid model. Also Section [3] states 
sufficient conditions for a fluid model solution to be unique, and for a fixed fluid model solution 
to be unique and asymptotically stable. Sections 2] and El discuss convergence of the scaled state 
descriptor and its stationary distribution to the fluid model and its fixed point, respectively. 
Sections El [7] and [8] contain the proofs of the statements from Sections El H] and El The Appendix 
proves auxiliary results. In the remainder of this section, we list the notation we use throughout 
the paper. 

Notation In order to introduce the notation, we use the signs =: and :=. 

The standard sets are denoted as follows: the reals R = (— oo, oo), the non-negative reals R+ = 
[0, oo), the positive reals (0, oo), the non-negative integers Z + = {0,1,2,...}, and the natural 
numbers N = {1, 2, . . .}. 

The signs A and V stand for minimum and maximum respectively. For x € R, x + := x V 0. 

The signs lim and lim denote the lower and upper limits of a sequence of numbers. 

The coordinates of a vector from a set S 1 are denoted by the same symbol as the vector with 
lower indices 1, ... ,1 added. If a vector has a superscript, tilde-sign, or overlining, they remain 
in its coordinates. For example x° € S d , x° = (x\, . . . ,x^). The space R 7 is endowed with the 
supremum norm ||x|| := maxi<j</ \xi\. Vector inequalities hold coordinate-wise. The coordinate- 
wise product of vectors of the same dimensionality I is x * y := (x\y\, . . . , xiyi). 

The signs =>, = and < s t stand for convergence in distribution, equality in distribution and 
stochastic order respectively. Recall that, for real-valued r.v.'s X and X' , X< st X' if ¥{X > x} < 
¥{X' > x} for all x € R. The notation 11(A), A € (0, oo), stands for the Poisson distribution 
with parameter A. 

For metric spaces S and S', denote by Cs->.s' the space of continuous functions /: S — > S'. By 
Ds + _>s denote the space of functions /: R + — > S that are right-continuous with left limits, and 
endow this space with the Skorokhod Ji-topology. 

The superscript —1 is only used to denote the inverse of a function. 

For a measure £ on R^_ and a £-integrable function /: — > R, define (/, £) := J r2 fd£. If 

£ = (£i, ...,£/) is a vector of such measures, (/, £) := ((/, • • • > ifiZi))- Let M be the space 
of finite non-negative Borel measures on R^_ endowed with the weak topology: £ fc — > £ in M as 
k — > oo if and only if (/, £ fc ) — > (f,Q for all continuous bounded function /: R^_ — > R. Weak 
convergence of elements of M is equivalent to convergence in the Prokhorov metric: for £, ip € M, 
define 

d(£ tp) := inf{e: £(B) < p{B e ) + e and cp{B) < £(B £ ) + e 
for all non-empty closed B C R^_}, 
where B £ = {x € R^_ : inf^^ ||x — y\\ < e}. 
For M 7 , define 

di(€,<p) ■= maxi<i</ d(£i, <pi). 
Equipped with the metric •), the space M 7 is separable and complete. 
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2 Stochastic model 



This section contains a detailed description of the model under consideration. In particular, it 
specifies the structure of the network, the policy it operates under and the stochastic dynamical 
assumptions. Also, a stochastic process is introduced that keeps track of the state of the network, 
see the state descriptor paragraph. 

Network structure Consider a network that consists of a finite number of links labeled by 
j = 1, . . . , J. Traffic offered to the network is represented by elastic flows coming from a finite 
number of classes labeled by i = All class i flows are transferred through a certain 

subset of links, we call it route i. Transfer of a flow starts immediately upon its arrival and is 
continuous with all links on the route of the flow being traversed simultaneously. Let A be the 
J x I incidence matrix, where Aji = 1 if route i contains link j and Aji = otherwise. 

Suppose that at a particular time t the population of the network is z S Z{, where Z{ stands for 
the number of flows on route i. All flows on route i are transferred at the same rate Xi(z) that 
is at most mi € (0, oo). If Z{ = 0, put Xi{z) := 0. We refer to Ai(z) := Xi(z)zi as the bandwidth 
allocated to route i. The sum of the bandwidths allocated to the routes that contain link j is the 
bandwidth allocated through link j and is at most Cj € (0, oo). We call Cj the capacity of link j. 
Hence, the vectors X(z) = (Ai(z), . . . , Xi(z)) and A{z) = (A\(z), . . . ,Aj(z)) must satisfy 

A(X(z) ■ z) = AA(z) < C, X(z)<m, A(z)<m*z, 

where C = (C±, . . . ,Cj) and m = (mj, . . . , m/) are the vectors of link capacities and rate 
constraints. 

Bandwidth-sharing policy At each point in time, the link capacities should be distributed 
among the routes in such a way that the network utility is maximized. Namely, to each flow 
on route i we assign a utility Vii(-) that is a function of the rate allocated to that flow. Assume 
that the functions Ui{-) are strictly increasing and concave in R + , and twice differentiable in 
(0, oo) with lima^oZ^(x) = oo. We also allow lim^o Ui{x) = — oo as, for example, in the case of 
a logarithmic function. Then, for z € R+, the vector X(z) of rates is the unique optimal solution 
to 

maximize / ZilAAXi) subject to A(X * z) < C, X < m, (1) 

where, by convention, x (— oo) := 0. Although the population vector has integer-valued coor- 
dinates, we assume that X(z) and A{z) := X(z) * z are defined via ([1]) in the entire orthant R^_ 
to accommodate fluid analogues of the population process later. 

The utility maximization procedure ([TJ implies that Xi{z) = Ai(z) = if z% = 0. The assumption 
lim^o^i^) = 00 guarantees non-idling, that is Xi(z), Aj(z) > if z\ > 0. Reed & Zwart [23] 
proved that the functions A(-) and A(-) are differentiable in any direction and, in particular, lo- 
cally Lipschitz continuous in the interior of R^. We also show continuity of A(-) on the boundary 
of R^ (see the Appendix). 

Lemma 1. The bandwidth allocation function A(-) is continuous in R^_. 

Stochastic assumptions All stochastic primitives introduced in this paragraph are defined 
on a common probability space (O, J 7 , P) with expectation operator E. 

Suppose at time zero there is an a.s. finite number of flows in the network, we call them initial 
flows. A random vector Z° £ R^_ represents the initial population, and Zf is the number of 
initial flows on route i. New flows arrive to the network according to a stochastic process £'(•) = 
(.Ei(-), ...,£/(•)) with sample paths in the Skorokhod space D R _^ R / . The coordinates of the 



4 



arrival process are independent counting processes. Recall that a counting process is a non- 
decreasing non-negative integer- valued process starting from zero. For t > 0, Ei(t) represents 
the number of flows that have arrived to route i during the time interval (0,i]. The kth such 
arrival occurs at time Ui k = inf{t > 0: Ei(t) > k}, it is called flow k on route i, k G N. 
Simultaneous arrivals are allowed. 

Flows abandon the network due to transfer completions or because they run out of patience, 
depending on what happens earlier for each particular flow. Flow sizes and patience times are 
drawn from sequences {(B^, Z^)}z eN , {(B ik , D ik )}ken, i = 1, of (0, oo) 2 -valued r.v.'s. For 
I = 1, . . . , Zf, B^ and Z)£ represent the residual size and residual patience time at time zero of 
initial flow I on route i. For k G N, Bik and Di k represent the initial size and initial patience time 
of flow k on route i, where "initial" means as upon arrival at time Un-. Let (B^, Di k ), k G N, be 
i.i.d. copies of a r.v. (Bi,Di) with distribution law 9i and a finite mean value (l//ij, Also 
introduce the notations pi := rjif/ii and Oi := rji/vi. Assume that the sequences {(B ik , D ik )} ke ^ 
are independent and do not depend on the arrival process E(-). For the moment, we do not 
make any specific assumptions about the sequences {(B^, D^) In- 



state descriptor We denote the population process by Z(-) = (Zi(-), . . . , Zj(-)), where Zi(t) 
is the number of flows on route i at time t. As can be seen from what follows, Z(-) is a random 
element of the Skorokhod space D R+ ^ R j . 

For i = 1, . . . , I, introduce operators Si : D K+ _ i , IR / x — >■ C^^^ defined by 

Si(z,s,t) — J Xi(z(u))du, 

For t > s > 0, Si(Z,s,t) is the cumulative bandwidth allocated per flow on route i during time 
interval [s,t]. The residual size and residual lead time at time t of initial flow / = 1,. . . ,Zf on 
route i are given by 

fl8(t):=(i?8-S(Z,(M))+ and £>g(t) := (Z?° - t) + , 

and those of flow = 1, . . . , Ei(t) on route i by 

flifctt) := (B ik - Si(Z, U lk ,t)) + and D ik (t) := (D ik - (t - U ik ))+. 

The state of the network at any time t is defined by the residual sizes and residual patience 
times of the flows present in the network. With each flow, we associate a dot in R+, whose 
coordinates are the residual size and residual patience time of the flow (see Fig. [1]) . As a flow is 
getting transferred, the corresponding dot moves toward the axis: to the left at the transfer rate 
(which is Xi(Z(t)) for a flow on route i) and downward at the constant rate of 1. As a dot hits 
the vertical axis, the corresponding flow leaves due to completion of its transfer. As a dot hits 
the horizontal axis, the corresponding flow leaves due to impatience. We combine these moving 
dots into the stochastic process Z(-) G D K+ _ i . M / with 

%i{t) ■= /Zi =1 S (B°(t),D^t)) + /Z k=1 S (B ik (i),D ik (t)y ( 2 ) 

where, for x G 5% G M is the Dirac measure at x if x\ A X2 > and zero measure otherwise 
(i.e. assigns a zero mass to any Borel subset of That is, Zi(t) is a counting measure on 
that puts a unit mass to each of the dots representing class i flows except those on the axis. 
The process Z(-) given by ([2]) we call the state descriptor. Note that the total mass of the state 
descriptor coincides with the network population, (1,Z(-)) = Z(-). 

When proving the results of the paper, we decompose the state descriptors into two parts keeping 
track of initial and newly arriving flows, respectively. That is, 

Z(-) = Z init (-) + Z ncw (-), 
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where 



Zi(0) 



and Zf™{t) :-- 



Em) , 
k=i (B ik (t),D ik (t)y 



We also define the corresponding total mass processes 

Z init (-) := (l,Z init (-)) and Z ncw (■) := {1, Z ncw (■)}. 
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Figure 1: The i-th coordinate -Zj(-) of the state descriptor puts a unit mass to 
the dots representing class i flows except those on the axis. 



3 Fluid model 

In this section we define and investigate a fluid model that is a deterministic analogue of the 
stochastic model described in the previous section. Later on the fluid model will be shown to 
arise as the limit of the stochastic model under a proper scaling. This convergence implies, in 
particular, existence of the fluid model. 

To define the fluid model we need data (77, 9, £°) € (0, 00) x x M J . The coordinates of r] 
play the role of arrival rates. Recall that 9i is the joint distribution of the generic size Bi and 
patience time Di of a newly arrived flow on route i. Finally, £° characterizes the initial state of 
the network. Also put z° := (1,C°) and, for all i, take a r.v. (Bf,D®) that is degenerate at (0,0) 
if zf = and has distribution Ci / z i otherwise. Then z° represents the initial population, and 
(B®, Df) the generic size and patience time of an initial flow on route i. 

Denote by C the collection of corner sets, 

C := {[x, 00) x [y, 00) : (x, y) G R 2 + }. 

Definition 1. A pair ((,z) € C K+ ^. M / x C K+ ^ K / is called a fluid model solution (FMS) for 
the data (n, 9, (°) if z(-) = (1, £(•)) and, for all i, t > and A € C, 

Ut)(A) =zM(B°,D°) € A + ($(*,(), t),i)} 

+ th [ ¥{(B i ,D i )eA+(S i (z,s,t),t-s)}ds. (3) 
J 
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In particular, for all i and t > 0, 

Zi (t) = Ci(t)O^) =*?P{B? > Si(z,0,t),D? > t} 

+ Vi I HBi>Si(z,s,t),Di>t-s}ds. (4) 
J o 

The function £(•) is called a measure-valued fluid model solution (MVFMS) and the function 
z(-) a numeric fluid model solution (NFMS) 

Equations ([3]) and @ have appealing physical interpretations. For example, simply means 
that a flow is still in the network at time t if its size and patience exceed, respectively, the 
amount of service it has received and the time that has passed since its arrival up to time t. 

Remark 1. By Dynkin's ir-\ theorem (see Section 2.3]), FMS's satisfy (J3J) with any Borel 
set A C R\. 

Remark 2. FMS's are invariant with respect to time shifts in the sense that, if (£, z){-) is an FMS, 
then, for any S > 0, (( s ,z s )(-) := ((,z)(- + 5) is an FMS for the data (rj, 9, £(£)). That is, for 
all i, i > <5 and Borel sets A C K+, 

Q(t)(A) = (i(5)(A + (Si(z, S, t),t - 5)) + m J ¥{(B U A) € A + (Si(*, s, i), t - s)}ds, (5a) 

z,(t) = &(6){[Si(z, 5, t), oo)x[t- S, 00)) + i lt j ¥{Bi > S t (z, s,t),Di>t- s}ds. (5b) 

Remark 3. The measure- valued and numeric components of an FMS uniquely define each other. 
In particular, uniqueness of an NFMS implies uniqueness of an MVFMS, and the other way 
around. 

Further we discuss sufficient conditions for an FMS to be unique and for an invariant (i.e. 
constant) FMS to be unique and asymptotically stable. To prove the stability result, we derive 
relations for asymptotic bounds for FMS's, which seems to be a novel approach since we have 
not seen analogous results in the related literature. We also give an example of multiple invariant 
FMS's. 

Uniqueness of an FMS The proof of the following theorem follows along the lines of the 
proofs of similar results [H Proposition 4.2] and Theorem 3.5], see Section 

Theorem 1. Suppose that either (i) 2? = for all i, or (ii) z° € (0, oo) 1 and the first projection 
°f C° * s Lipschitz continuous, i.e. there exists a constant L € (0, 00) such that for all i, x < x' 
and y, 

$([x,x'] x [y,oo)) <L{x'-x). 
Then an FMS for the data (rj,6,C°) ^ unique. 

Uniqueness of an invariant FMS Let (£, z) be an invariant FMS. By Lemma[3]in Section H3 
all of the coordinates of z are positive, and the fluid model equations (|3D and flU for (£, z) look 
as follows: for all i, Borel subsets A C M. 2 ^ and t > 0, 

Ci{A) = Ci(A+(X i (z)t,t))+r Ji f 6i(A + (Xi(z)8,s))ds, (6) 

J 

Zi = <i([\i(z)t, 00) x [t, 00)) + Vi / F{Bi > \i(z)s, A > s}ds. (7) 

J 
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Letting t — > oo in (JSJ) and ([7]), we obtain the equations 

/•oo 

Ci(A) = Vi / 0i(A + (Xi(z)s, s))ds, (8) 
Jo 

z i = r H E(B i /X i {z)AD i ), (9) 
which are actually equivalent to (|6|) and ([7]). 

Thus, we have the closed-form equation ([9]) for the numeric components of invariant FMS's, and 
the corresponding measure- valued components are defined by ([8]). In particular, uniqueness of 
an invariant FMS is equivalent to uniqueness of a solution to ([9]). 

Definition 2. For an M- valued r.v. X, denote by inf X the left most point of its support. Recall 
that the support of X is the minimal (in the sense of inclusion) closed interval 5 such that 
¥{X € S} = 1. 

Theorem 2. Let > inf(Dj/ ' Bi) for alii. Then there exists a unique invariant FMS z*), 
and the bandwidth allocation vector A(z*) is the unique solution to the optimization problem 

maximize ZiGAKi) subject to AA<C, A, < 7]iE(Bi A rriiDi) for alii, (10) 

with strictly concave functions Gi(-). 

To prove the theorem, we use the approach of Borst et al. [U Lemma 5.2], which is as follows. For 
now, let z* be any invariant NFMS As we plug the fixed point equation ([9]) into the optimization 
problem (fTJ) for the rate vector \(z*), the problem (|1U|) follows, which is strictly concave and 
does not depend on z*. Hence, A(z*) is the same for all invariant points z*. The assumptions of 
the theorem guarantee invertibility of the function A(-), and then uniqueness of A(z^) implies 
uniqueness of z*. 

Note that this method not only proves uniqueness of an invariant FMS, but also suggests a way 
to compute it (a strictly concave optimization problem can be solved with any desired accuracy 
in a polynomial time). 

Asymptotic bounds for FMS's Here we derive asymptotic bounds for NFMS's that, for a 
wide class of bandwidth-sharing networks, imply convergence to the invariant NFMS provided 
it is unique. 

Theorem 3. There exist constants l,u G (0,00)^ such that, for any NFMS z(-), 

< k < lim^^Ziti) < lim^oo^) < m for all i. 
These constants satisfy the relations 

l i = r H E(B i /R i {l,u)AD i ), Ui = rnE(Bi/ri(l, u) A A) for alii, (11) 
where the functions r(-, •) and R(-,-) are defined by 

ri{x,x) := mi x < z < x i \i(z), Ri(x,x') := sup x<z<x i Xi(z) for all i and x < x . 

Remark 4. There could be more than one pair (/, u) solving (|lip . The asymptotic bounds / and u 
for NFMS's given by Theorem [3] form one of such pairs. 

We now proceed with the proof of Theorem [3l 
Proof. Note that if 

< k < \im t ^^Zj(t) < limt^ooZiit) < Ui for all i, (12) 
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then also 

rjiEiBi/Rid u) A A) < Hm^oo^W < Enw^i) < r]iH B i/n{l u) A A) for all i. (13) 
Indeed, by (fT2|) . for any e € (0, mini<j</ Zj), there exists a i e such that 

li — e < Zi(t) < Ui + £ for all i and t > t e . 

Introduce the vectors 

I — £ := (h — £, . . . , li — e), it + e := (£ti + e, . . . , itj + e). 

Then 

Ti(l — £,u + e)(t — s) < Si(z, s, t) < Ri(l — £,u + e)(t — s) for t > s > t £ , 
which, when plugged into the shifted fluid model equation (|5b|) . implies that 

Zi(t)>TH t F{B i >R i (l-£,u + £)(t-s),D i >(t-s)}ds, 

Jt e 

Zi(t) < d(t e )([Si(z,t £ ,t),oo) x [t-t e ,oo)) 

+ rji f F{B i >r i (i-£,u + £)(t-s),D i >(t-s)}ds for t > t e , 

Jte 

where C(-) is the corresponding MVFMS. Taking t — > oo in the last two inequalities, we obtain 

rtiE(Bi/Ri(l-£,u + e) A A) < lim^ ~z«ffl < lim^oo^^) < i) i E(B ! /r ! ([ - e, it + e) A A), 
and then (|13p follows as e — )• 0. 

Now we will iterate f)12[> — f)13|) . The rate constraints plugged into (j3|) imply the initial bounds 

< if := rtiE(Bi/mi A A) < lim f ^Zi(t) < lim^^^^) < ^EA =: t£ for all i, 
and then (|12|) — (jl3 [) yield the recursive bounds 

l\ := nt E(B t /R t (l k -\u k - 1 ) A A) < lim^^i), 

, , , , , (14) 

< := r] i E(B i /r i (l k - 1 ,u k ~ 1 ) A A) > lim^ooZ^t) for all k G N and i. 

The sequence {/ fc }fceN is non-decreasing and bounded from above by vP. The sequence {u k }k^n 
is non-increasing and bounded from below by 1°. Hence, there exist the limits lim/ fc =: / and 
limu k =: u. In let — >■ oo, then (fTT]) follows. 

Note finally that the recursive bounds {^jfceN and {u fc }fc g N as well as their limits I and u do 
not depend on a particular NFMS. □ 

Asymptotic stability of an invariant fluid model solution It is realistic to assume that 
transfer rates in a bandwidth-sharing network decrease as its population grows. 

Definition 3. If z' > z € (0, oo) implies X(z') < X(z), the network is called monotone. 

For monotone networks, the system of equations (|lip decomposes into the fixed point equa- 
tion (J9j) for the lower bound I and for the upper bound u, implying the following result. 

Corollary 1. Suppose that the network is monotone and has a unique invariant FMS (£*,z*). 
Then any FMS (C, z)(t) ->• ((*, z*) as t ->• oo. 

(It is easy to see how Theorem [3] implies z(t) — > z* in Corollary [TJ In the Appendix we also show 
that z(t) ->• z* implies £(i) -> C*0 

Example: single link & multiple classes The sufficient condition for uniqueness of an 
invariant FMS given by Theorem[2]is sometimes also necessary. Consider, for example, processor 
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sharing in critical load, that is J = I = 1 and (omitting the link and class indices) p = C. In 
this case, the fixed point equation ([9]) looks like 



which, for z such that C/z < m and Bz/C < D a.s., reduces to 

z = 7]E(Bz/C A D) = rjEBz/C = pz/C = z. 

I.e. any z € [C/m, Cinf D/B] is an invariant NFMS. In particular, if m.iD/B > 1/m, which 
violates the assumption of Theorem [21 then there is a continuum of invariant FMS's. 

For a single link critically loaded by multiple classes of flows, we have an analogous result, which 
is more complicated to derive and therefore the proof is postponed to Section [6J 

Theorem 4. Assume that J = 1 (in what follows we omit the link index), and that the utility 
functions are lAi{x) = Kjlogx. If Yli=i Vi^(Bi/ m i A Ai) C> then there is a unique invariant 
FMS. Otherwise there might be a continuum if invariant FMS's. 

4 Sequence of stochastic models and fluid limit theorem 

In this section we study the asymptotic behavior of the stochastic network described in Section [2] 
as its global parameters — capacities and arrival rates — grow large, while the characteristics of 
an individual flow remain of a fixed order. We refer to this scaling as the large capacity regime. 

Large capacity scaling To a sequence 7Z of positive numbers increasing to oo, we associate 
a sequence of stochastic models as defined in Section [21 We mark all parameters associated with 
the r-th model with a superscript r and assume the following: 

(A.l) network structure, rate constraints and utility function are the same in all models: 
A r = A, m r = m and U[(-) = Ui(-) for all i; 

(A. 2) link capacities grow linearly in r: C r = rC; 

(A. 3) arrival rates grow linearly in r: 'E' (■) := E r (-)/r =>■ n{-) as r — > oo, where n(t) := trj and 



(A. 4) flow sizes and patience times remain of a fixed order: for all i, (Bl,D\) (Bi,Di), 
where (Bi,Di) are (0, oo) / -valued r.v.'s with distributions 9{ and finite mean values 
(1/fH, l/vi), and also 1/z-f ) -> (1/fH, 1M); 

(A. 5) the scaled initial configuration converges in distribution to a random vector of finite 
measures: if'(O) := Z r (0)/r (°; 

(A. 6) the projections x K + ) and Cf(^+ x ") are a - s - f ree °f atoms for all i. 

Fluid limit theorem In the large capacity regime, the stochastic model defined is Section [2] 
converges to the fluid model defined in Section [3J More precisely, introduce the scaled state 
descriptors and population processes 




V G (0,oo) 



r (•) := Z r (-)/r, Z r (-):=(l,Z\-)) = Z r (.)/r. 



Also introduce the scaled versions of the two components of the sate descriptor: 



Z r ' ,mit (-) := Z r ' init (-)/r, T' (■) := (1, Z r ' imt (-)> = Z r ' init (-)/r, 

Z r ' ncw (-) := Z r ' ncw (-)/r, Z r ' ncw (.) := (l,Z~ r ' ncw (■)) = Z r ' ncw (-)/r. 
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Remark 5. Let A(-) be the rate allocation function in the unsealed network, then 

\ r (z) = X(z/r), 

S^(Z r ,s,t) := f \\{Z r (u))du = I \i(Z r (u))du=: Si(Z r ,s,t). 

J s J s 

We now provide the definition of a fluid limit followed by the main result of this section. 

Definition 4. We refer to weak limits along convergent subsequences {(Z q , Z 9 )(-)} (?e g, Q C 7£, 
as fluid limits. 

Theorem 5. Under the assumptions (A.1)-(A.6), the sequence {(Z r , Z r )(-)} r£ -ji is tight in 
D K+ _ >M / x D R+ _ S>K ^, and all fluid limits are a.s. FMS's for the data (j],9,(°). In particular, 

if there is a unique FMS (Z, Z)(-) for the data (rj, 9, (°), then [Z\ Z r )(-) => (Z, Z){-) as r — > oo. 

The proof follows in Section [JJ To show tightness we adjust the techniques of [13] to the two- 
dimensional case, since in [13] flows are patient and state descriptors are vectors of measures on 
R+. The proof of convergence to FMS's follows the lines of that in [llj . It uses the boundedness 
of fluid limits away from zero, and the key difference is that in [llj this property is guaranteed by 
the overload regime, while in our model it holds in any load regime due to the rate constraints. 

5 Convergence of stationary distributions 

Assume that, in the stochastic model defined in Section [21 the arrival processes are Poisson of 
rates r]i,...,i]i. Then there exists a unique stationary (and also limiting as t — > oo) distribution 
of the state descriptor Z(-). Indeed, without loss of generality, there are i.i.d. r.v.'s {Dik}keN,i<i<i 
distributed as maxi<j</ Di and such that a.s. < for all k and i. Then the total population 
Y^l=i Z-i(-) of the network is a.s. and within the whole time horizon R + bounded from above 
by the length of the M/G /oo queue with the following parameters. At time t = 0, there are 
X)j=i Zi(0) customers in the queue whose service times are patience times of initial flows in 
the network. The input process for the queue is the composition of those for the network, and 
hence is Poisson of rate Yli=i Vi- Service times of new customers in the queue are drawn from 
the sequence {Dik} fceN,i<«</ of upper bounds for patience times of new flows in the network. 
As any other M/G/oo queue, the defined queue is regenerative. The instants when a customer 
enters the empty queue form an embedded renewal process whose cycle length is non-lattice and 
has a finite mean value exp(V\ =1 rji ^Dn)/Yli=i Vi- With respect to this renewal process , the 
state descriptor Z(-) is also regenerative. Then, by [H Chapter V.I, Theorem 1.2], there exists 
a limiting distribution for Z(-). 

Now consider a sequence of stochastic models as defined in Section [2] that satisfies the assump- 
tions (A.l), (A.2), (A. 4) (see Sectiong]) and 

(A'. 3) the input processes E\ (•),..., Ej(-) are independent Poisson processes of rates 77 J , . . . , 77J , 
and rf jr — >■ r\ G (0, oo)^ as r — > 00. 

Let y r have the stationary distribution of Z r {-) and put Y r := (1,3^). Introduce also the scaled 
versions 

y:=y r /r, T :={l,y r )=Y r /r. 
We now have the following result. 

Theorem 6. Under the assumptions (A.l), (A.2), (A'. 3) and (A. 4), the sequence {(y r ,Y r )} re fi 
is tight, and any weak limit point (3^,^) is a weak invariant FMS, i.e. there exists a stationary 
FMS (Z,Z)(t) = (y, Y), t > 0. In particular, by Corollary [7J if the network is monotone and 
has a unique invariant FMS (£*,£*), then (y ,Y ) (£*, z*) as r — > 00. 
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The general strategy of the proof is adopted from \17\ Theorem 3.3]: we check that any con- 
vergent subsequence of initial conditions {Z q (0) = y 9 } q eQ, Q C TZ, y q => y, satisfies the 
assumptions of the fluid limit theorem (we only need to check (A.6)). Then the correspond- 
ing subsequence {Z {-)} q eQ of the scaled state descriptors converges to an MVFMS that is 
stationary (i.e. 3* is a weak invariant MVFMS) since all Z q {-) are stationary. 

The techniques we use to implement this strategy are different from the techniques of [17] . 
though. Our key instruments for establishing tightness are M/G/oo bounds, see Section [HJ 
Below we present an elegant proof of (A.6) for weak limit points of {3^ }reK- 

Lemma 2. Any weak limit point y of {y r }ren has both projections y(- x K + ) and y(M.+ x •) 
a.s. free of atoms. 

Proof. The key idea is the following. Consider the network in its stationary regime. Then, on 
one hand, it always has the same distribution, and on the other hand, all initial flows are gone 
at some point, and newly arriving flows do not accumulate along horizontal and vertical lines. 

Let 3^ be the weak limit along a subsequence {y q }qeQ, and run the q-th network starting from 

Z q (0) = y q . By pTJ Lemma 6.2], it suffices to show that, for any 5 > and e > 0, there exists 
an a > such that 

Hm^PHsup,^ ||3^ +a )ll V ||3^ +a )|| < 5} > 1 - e, (15) 
where H b a : = R + x [a, b] and V a fe := [a, b] x R + for all b > a > 0. 

First we estimate the time when there is only a few (when scaled) initial flows left. The initial 
flows whose initial patience times are less than t are already gone at time t. Then Lemma [TBI (see 
Section [8|) implies that (recall that 11(A) stands for the Poisson distribution with parameter A) 

Zf ni \t) <Z q (0)(R + x [t,oo))±y q (R + x [i,oc)) 

1 f'OO f'QO 

< st / P9{U t ? > y}dy) m / P{A > y}dy. 

Q Jt Jt 

Take T such that maxi<j</ rji f°° P{-Dj > y}dy < 5/2, then 

lim q ^ 00 F r {\\Z q ' ini \T)\\<5/2} = l. 

Now, in Lemma [11] (see Section [7]), we prove that newly arriving customers do not accumulate 
in thin horizontal and vertical strips, i.e. there exists an a > such that 

hm^ooP^sup^r] sup xm+ \\Z q > ncw (t)(H* +a )\\ V ||^' new (t)(^ +a )|| < 5/2} >l-s. 
Finally, because of stationarity of y q , 

v q {su Pxm+ \\y q (H* +a )\\ v \\y q (v* +a )\\ < 5} 

=F«{sup xeR+ \\Z q {T){H^)\\ V \\Z q (T)(V* +a )\\ < 5} 

>F0{\\Z q '' mi \T)\\ < 5/2,sup xm+ \\Z q ' mi \T)(H^ a )\\ V p 9,init (r)(^+°)|| < 5/2}, 
which implies (|15p by the choice of T and a. □ 

6 Proof of fluid model properties 

Here we prove the results of Section [3j 
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6.1 Proof of Theorem [T] 



In the proof of Theorem [TJ we exploit boundedness of NFMS's away from zero and in the norm 
(see Lemma [31), and Lipschitz continuity of MVFMS's in the first coordinate (see Lemma [5]). We 
also use the auxiliary Lemma [TJ 

Lemma 3. Let z(-) be an NFMS Then sup t>0 ||z(i)|| < 11-2(0)11 + ll°"ll < 00 and, for any 5 > 0, 
inf t >5 mini<j</ Zi(t) > 0. In particular, if ZiiO) > 0, then 'wIt>oZi(t) > 0. 

Proof. By the rate constraints, Si(z, s, t) < mj(t — s) for all s < t, which, when plugged into the 
fluid model equation (|1J), implies the following lower bound: z.i(t) > rji f F{Bi/rrii AD; > s}ds. 
Since fi(s) := ¥{Bi/rrii A Di > s} t P{-Bi/mi A Di > 0} = 1 as s J, 0, in a small enough interval 
[0,e], fi(-) > 1/2. Then, for t > 5, Zi{t) > m J Q Me fi(s)ds > (5 A e)/2. The upper bound follows 
from (gD directly: Zi{t) < ^(0) + 7ft J * P{A > s}ds t ^(0) + cr, as i 1 °o. □ 

Lemma 4. For an M-valued r.v. £ and x < x', J R P{u + x < £ < u + x'}du < x' — x. 
See the Appendix for the proof. 

Lemma 5. Under assumption (ii) of Theorem [JJ any MVFMS Q{-) at any time t > has 
a Lipschitz continuous first projection, i.e. there exists a constant L(£,t) € (0, oo) such that for 
all i, x < x' and y, 

Ci(t){[x,x'} x [y,oo)) <L((,t)(x' -x). 

Proof. For an FMS (£, z)(-), for all i, i > 0, x < x' and y, 

d(t)([x,x'} x [y,oo)) < fi(x,x',y) + 7]igi(x,x',y), 

where 

fi{x,x',y) := C^([x + S l (z,0,t),x' + S l (z,0,t)] x [j/ + t,oo)), 

gi(x,x',y) := ¥{x + Si{z, s,t) < B { < x' + Si(z, s,t)}ds. (16) 

By Lipschitz continuity of the initial condition, fi(x,x',y) < L{x' — x). In (|16[) . change the 
variable of integration for v = V(s) := Si(z,s,t). Then 

rSi(z,0,t) 

gi (x, x', y) = F{x + v < B % < x' + v} / h^V' 1 {v)))dv < Af (£ - x), 

JO 

where M(£f) := sup sg r ot i maxi<j</ 1/Xi(z(s)). By Lemma[3l the functions 1/A;(z(-)) are con- 
tinuous in [0,t]. Hence M(£,t) is finite and the first projection of £(t) is Lipschitz continuous 
with the constant L((, t) := L + \\r)\\M((, t). □ 

Now we are in a position to prove Theorem [TJ 

Proof of Theorem^ Let (C 1 ,^ 1 )(-) and (( 2 ,z 2 ){-) be two FMS's for the data (n,6»,C°). 

(i) We show that the two FMS's coincide in an interval [0,5]. We check that z 1 (5) = z 2 (5) £ 
(0, oo) 1 and that the first projection of C 1 ^) = C 2 (^) is Lipschitz continuous. Then, by Remark[5] 
and the second part of the theorem, the two FMS's coincide everywhere. 

Note that, for a vector z G (0, oo) 1 of a small enough norm, \{z) = mi for all i. Lemma [3] and 
the fluid model equation (jTJ) imply that < zj(t), zf(t) < Tfti for all i and t > 0. Then, for all i 
and s, t 6 [0,5], where 5 is small enough, 

Si(z l , s, t) = Si(z 2 ,s, t) = rm(t - s). (17) 
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Plugging ()17p into (J1J), we obtain, for t £ [0, 5] and all i, 

4(i) = zf{t) =rji [ F{B t / mi A A > 

JO 

By Remark [3l C 1 (') an d C 2 (0 coincide in [0,5], too. Lipschitz continuity of the first projection 
of C 1 ^) = C 2 (<5) follows as we plug (fTT|) into the fluid model equation ([3]) (recall that it is valid 
for all Borel sets): for all i, x < x' and y, 

d(6)([x,x'] x [y,oo)) = rji / ¥{x + rriiS < Bi < x' + rriiS, A > U + s}ds 

Jo 

,6 

< rji / P{x/mj + s < Bi/rrii < x /mi + s}<is 

< rn(x' - x)/nii, j = 1,2, 
where the last inequality holds by Lemma HI 

(ii) Suppose that the two FMS's are different, that is i* := inf{i > 0: z it) ^ z 2 (t)} < oo. 

Without loss of generality we may assume that t* = 0. Indeed, otherwise we can consider the 
time-shifted FMS's (£ J , z 3 )(t* + •), j = 1,2. By Lemmas [3] and O they start from = 
z 2 (t*) € (0,oo) 7 and C X (**) = C 2 (^*) with a Lipschitz continuous first projection. 

By LemmaEl the two NFMS never leave a compact set [5, A] 1 C (0, oo) . Since the rate functions 
Xi(z) are Lipschitz continuous in such sets, there exists a constant K £ (0,oo) such that, for 
all i and s <t, 

\Si(z\s,t) -5i(z 2 ,s,t)| < Ktsup se[0>t] Hz 1 ^) -z 2 (s)|| =: iTte(t). 
Then, by Lipschitz continuity of the initial condition, we have, for all i and t > 0, 

|4 (t) - z?(t)| < LKte{t)+r]i [ P{5 i (z 1 ,s,t) - Kte(t) < Bf < Siiz 1 , s,t) + Kte{t)}ds. 

Jo 

In the last equation, change the variable of integration for v = Si(z 1 ,s,t) (c/. the proof of 
Lemma|5| and put M = sup^r^U max i<i</ l/\i(z). Then 

\z}{t) - zf{t)\ < LKte{t) + r]iM2Kte{t) for all i and e(i) < (L + 2||r/||M)iiCte(i), 

which implies that e(t) = for small enough t, and we arrive at a contradiction with = 0. □ 

6.2 Proof of Theorem H 

Multiplying the coordinates of ([9]) by the corresponding rates Xi(z), we obtain the equivalent 
fixed point equation 

Hz) = gi (Xi(z)) foralH, (18) 

where 

5i(x) := r]iE(Bi A xA), x > 0. 

We first discuss some properties of the functions gi{-) in the auxiliary Lemmas [6] and [71 and then 
proceed with the proof of the theorem. 

Lemma 6. The function gi{-) is continuous. Also gi{-) is strictly increasing in [0, «j) and con- 
stant in \ [0,ai), where 

Oi := inf{x: g^x) = pi} > 0, 
and infimum over the empty set is defined to be oo. 

Proof. Continuity of <?j(-) follows by the dominated convergence theorem. 



14 



The situation ctj = is not possible since in that case gi(x) = p% for all x > by the definition 
of on. But gi{-) is continuous and g%{x) — > gi(0) = as x — > 0. 

If aj < oo, then, again by the definition of a, and continuity of <?i(-), we have gi{x) = pi for all 

X > OLi. 

It is left to check that gj(-) is strictly increasing in [0, aj). Assume that < x < y < oti, but 
Si 0*0 = <7i(y)- Then 

o = gi{y)/vi - gi{x)hi =EBil {B .< xD .j +E5 i i {:rDi<i j i < yjDi} + EyA^s^^} 

- EBjl^.^^o.} - ExDiI {xDi<BiSyD . } - ExDiI {Bi>yDi} ) 

=E (J3» - »A)I{xD i <B i <y£» i } +(?/ - X)E DiI {B . >yD .y, 
s v ' v v ' 

= :X =:Y 

where the r.v.'s X and Y are non- negative, so they must a.s. equal zero. In particular, we have 
Bi < yDi and gi{y) = Pi, which contradicts the definition of «j since y < on. □ 

The stabilization points on of the functions gi(-) are related with the r.v.'s (Bi,Di) in the 
following way. 

Lemma 7. If on < oo, then vniDi/Bi = l/ctj. = oo, then miDi/Bi = 0. 

Proof. First assume a, < oo. Rewrite the relation gi{x) = pi as E5j(l — (1 A xDi/Bi)) = 0, 
which, for x > 0, is equivalent to Di/Bi > 1/x a.s. Hence Oj = inf{x > 0: Di/Bi > 1/x a.s.} 
and 1/ctj = sup{y > 0: Di/Bi > y a.s.}. In the right-hand side of the latter equation we see the 
definition of inf Di/Bi. 

Now consider the case on = oo. Assume that inf Di/Bi = y > 0, then Di/y > Bi a.s. and 
gi(l/y) = pi. On the other hand, since on = oo, there is no x > such that gi(x) = pi. Hence 
y = 0. □ 

Having established the above properties of ^ (-)'s, we now can prove Theorem [2] by means of a 
technique developed by Borst et al. [H Lemma 5.2]. 

Proof of Theorem^ First show uniqueness. Let z* 6 (0, oo)^ be an invariant NFMS. Recall that 
X(z*) is the unique optimal solution for the concave optimization problem ([1]). The necessary 
and sufficient conditions for that are given by the Karush-Kuhn- Tucker (KKT) theorem (see 
e.g. [31 Theorem 3.1]): there exist p G and q € such that 

U'i(h{z*)) =Y, J j=1 A JiPj + <H for alii, (19a) 

p j (^ I . =i A ji k i {z*)-C j )=U for all j, (19b) 

qi(Xi(z*) — rrii) = for all i. (19c) 

By the assumptions of the theorem and Lemmas [6] and [TJ the functions gi(-) are strictly in- 
creasing in the intervals [0,mj], which implies two things (see also Fig. [2]). First, the fixed point 
equation (fTHj) can be rewritten as \{z*) = g J ~ 1 (Aj(z*)) for all i, and we plug that into (|19aD . 
Second, the second multiplier in (|19c|) is zero if and only if gi(Xi(z*)) = (jj(mj), and that, by (fTHj) . 
is equivalent to An(z*) = gi(mi). Hence, A(z*) satisfies 

Ui{gr\K t (z*)) = Y? A jiPj + q* for all i, (20a) 

PjC£! i=1 A iMz*)-Cj) = for all j, (20b) 
q i {k i {z*)-g i (m i )) = U for alH. (20c) 
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Now note that the last three equations form the KKT conditions for another optimization 
problem. Indeed, take functions g%{-) that are continuous and strictly increasing in R + and 
coincide with gi(-) in [0, m^]. Also take functions Gj(-) such that G[{-) = U^g^ 1 ^)) in (0, oo). 
Then (|20p gives necessary and sufficient conditions for A(z*) to 

maximize ZiGAkj) subject to AA<C, A» < gdrrii) for all i. (21) 

Since the functions Ui{-) are strictly concave, their derivatives U[{-) are strictly decreasing. Then 
G[{-) are strictly decreasing and, equivalently, Gj(-) are strictly concave, which implies that 
A(z*) is actually the unique solution to (|21|) . Then the invariant NFMS z* must be unique, too, 
z* = Ki{z*)/gr l (hi{z*)) for aU i 

The existence result follows similarly. There exists a unique optimal solution A* to (|2ip and it 
satisfies the KKT conditions ([2D]). Put A* = ^(A*) for all i. Then A* and A* satisfy the KKT 
conditions ([12]), i.e., for the vector z* with z* := A|/A*, we have A* = X(z*) and A* = A(z*). 
Plugging the last two relations into the definition of A*, we get the fixed point equation. □ 




Figure 2: Graph of the function gi(-) in the two possible cases: 
when rrii < a, (left) and when m« > (right); z is a invariant NFMS. 



6.3 Proof of Theorem H 

The fixed point equation (|18j) and the monotonicity of the functions gi{-) imply that the band- 
width class i gets in an equilibrium is at most gi(rrii). Therefore, we refer to the scenarios 
£i=i 9i( m i) < Yli=i9i( m i) = C and ^i=i5«( m «) > C as underloaded, critically loaded and 
overloaded, respectively. Below we calculate the invariant NFMS's in the three cases. 

Summing up Q, the KKT conditions for ([T]) and the capacity and rate constraints, a z € (0, oo) 1 
is an invariant NFMS if and only if there exist p G M+ and q G such that (we omit the 
argument of the rates \%{z) and bandwidth allocations Aj(z)) 



A l =g i (X i ) for alii, (22a) 

= p + q% for all i, (22b) 

p(V J A - C) = (22c) 
* — »i=i 

qi{Xi - mi )=0 for all i, (22d) 

V 7 A, < C, (22e) 
* — 'i=i 

A» < mj for all i. (22f) 



Underload In this case, there is no interaction between the classes, they do not compete but 
all get the maximum rate allowed. Indeed, ()22cp and (|22bp imply that p = and all gj > 0. 
Then, by ()22dp and ()22ap . all Aj = rrij and all A» = gi(nii). Hence, there is a unique invariant 
NFMS given by 

Zi = gi(mi)/mi for all i. 
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Critical load First note that 

Aj = gi{mrii) for all i. (23) 

Indeed, there are two possibilities: either p = 'W^' all > 'W"'' all Aj = m.j (|2"3"1) . or p > 

=P X^f = i Aj = C => (|23|) . where the last implication is due to Aj < gj(m-j) and 9i( m i) = C- 
Recall from Lemma [7] that 

Qi := inf{x: g^x) = pi} = 1/inf (A/A)- 

By (|2"3"j) . the relations (|22a|) and ()22f|) are equivalent to m. L A ctj < Aj < mj ( see Fig. [2]). Hence, 
([22]) reduces to 

Ki/\i=p + qi, (24a) 
g i (A i -m i ) = J (24b) 
m; A Qj < Aj < mj. (24c) 

Let 

Zcrit := {i: mj < aj}. 

For i € Xcrit, by (|24c|) . we have Aj = mj and m = <7i(mj) /mj. Then (|24bp is satisfied, and (|24a|) 
implies that p < Kj/mj. 

Now divide {1, ...,/} \ X cr ;t into two subsets X\ n X2 — 0. For i € Ii, put Aj — mj, then (as 
for i € 2crit) 2j = <7j(mj)/mj, ()24bp is satisfied, and (|24ap implies that p < ftj/mj. For i G I2; 
assume aj < Aj < mj. Then <?j = by (|24bp . ftj/Aj = p by (|24ap . and Ki/rrii < p < Hi/mi. Also 
z = gi(mi)p/Ki. 

Summing up everything said above, the set of invariant NFMS's is given by 

S z : = M {z: Zi= gi(mi)/mi for i G X and Zj = g^m^p/Ki for i ^ X, 

where p € (maxj^ Kj/mj, min ieI Kj/mj A mim.^ Hi/an]}. 

Equivalent descriptions of S z are 

S z = {z: Z{ = gi(mi)/mi if p < Kj/mj and Zj = gi(rrn)p/Ki if p > Kj/mj, p 6 
:= (0, minj e i crit Kj/mj A mhij^^ Ki/oci] = (0, mini<j</ m/ (jm A a;)], 

and 

S z = {z: Zi = gi(mi)/(mi A Kix), x £ S x }, 
S x := [maxi<j</(?nj Aaj)/Kj,oo). 

We now apply the last formula in a couple of simple examples. 

Example 1. If mj < aj for all i, then S x = [maxi<j</ mj/ftj, 00), and K{X > mj for all x € S x and 
all i. Hence, there is a unique invariant NFMS given by Z{ = giijn^jmi for all i, which agrees 
with Theorem ([2]). 

Example 2. If m\ > a\, n%i < Oj for i 7^ 1 and > maxj^x mj/^j, then, for any Ai € mi], 

z = (5 , i( m i)/^i)fi , 2(^-2)/' 7 i2; • • • ^gi{jni)jm{) is an invariant NFMS. 



Overload In this situation, by the capacity constraint (|22e|) . at least one class of flows does 
not receive the maximum service, i.e. at least one Aj < <7i(mj). We first find out which classes 
get the maximum service and which do not, and then calculate the unique invariant NFMS. 

Who gets the maximum service. Since at least one Aj < ft(mj), at least one Aj < mj A aj (see 
Fig. [2]). Then (j22d|) . (f22"b"|) and (l22c|) imply that at least one % = 0, p > and Y!i=i A i = C. 
At this point, we can equivalently rewrite (|22p as follows: there exist x > and e S such 
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that (the functions <ji{-) are introduced below) 

K = 9i{\i) Ai = (25a) 
, ^(Ai) = C & ft(Aj) = C (25b) 

i — *l=\ £ '4 = 1 

\i = Ki(x-£i), (25c) 
£i (Ai - mi) = 0, (25d) 
Ai < mi. (25e) 

For all % and x > 0, put 

gi(x) := gi(mi A x). 

By the rate constraints (|25ej) . in ([25]) . we can equivalently replace gi(-) by gi(-). 
If Aj < m;, then, by (j25dp and (|25cp . £j = and Aj = KiX, and hence 

3i(Ai) = 9i{Kix). (26) 
If Aj = mi, then, by (|25c[) . KiX > mi and <?j(«;jx) = ft(mj), and, again, ([26]) holds. 
Plugging ([25]) into (|25b|) . we get 

V 7 = C. (27) 

* — ^j=i 

The function g(x) := Yli=i9i( K i x ) ^ s continuous everywhere, strictly increasing in the interval 

< x < maxi<j</(?nj A on)/ Ki =: xq 

and constant for x > xq, and also g(0) = and g(xo) = Yll=i9i{ m i) > C) which implies that 
there exists a unique x solving ([27]) and x G (0,xo). 

By ()25ap and ([26]) . Aj = cn(Kix). Then (see Fig. [3]) Aj = gi(m,i) if (m; A a,i)/ni < x and 
Aj < gi(mi) if (mj A aj)/Kj > x. Hence, the set of classes that get the maximum service is 

lover := {i- (m» A Oi)/«i < x}. (28) 

Invariant NFMS. For i ^ X over , Aj = gi(Xi) < <7i(m.j), which implies that (see Fig. [2J Aj < 
mi A ctj. Then, by ()25dp and (|25cp . e$ = and Aj = KjX (meeting the rate constraint (|25ep ). and 
Zi = Ai/Xi = gi(Kix)/(Kix). 

For i € X over , consider the two possible cases: KjX < m; and KjX > mj. If KjX < mj, then, by (|25cp 
and (j25dp . Aj < KjX < mj and £j = 0, and, again by (|25cp . Aj = re,x. If k^x > mj, then Aj = mi 

because otherwise we would arrive at a contradiction: Aj < mi e% = Aj = fqx > mj. 
Hence, for i 6 X over , Aj = mj A KjX and Zj = Aj/Aj = gi(mi) / (mi A Kjx). 

Summing up, the unique invariant NFMS is given by 

Zi = gi(mi)/(mi A Kjx) for i G X over and Zi = ^i(«tx)/(/Cfx) for i ^ Z over , 
where x is the unique solution to ([27]) and 2o Ver is defined by 



7 Proof of Theorem [5] 

To prove C-tightness (that is, tightness with all weak limits being a.s. continuous) of {Z r (-)} r £TZ, 
we check standard conditions (see e.g. [10]) of compact containment (Section 17. 2p and oscillation 
control (Section l7.4p . In Section lT^Bl we check that fluid limits satisfy the fluid model equation ([3]). 

To establish these main steps of the proof, we develop a number of auxiliary results. Section I7TT1 
contains a law of large numbers result for the load process. Section [7.31 proves that, for large r, 
Z (•) puts arbitrarily small mass to thin horizontal and vertical strips, which in particular implies 
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9j(mj) 
g l {m i )=g l {n t x) 
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1 * 





(a, A m,i)/Ki x (a.j A rrij)/ Kj 

Figure 3: Graphs of the functions &(■); x is the unique solution to ([2 



that fluid limits have both projections free of atoms. In Section T7.5| fluid limits are shown to be 
coordinate-wise bounded away from zero outside t = 0. 

7.1 Load process 

Introduce the measure valued load processes and their scaled versions: for all r, i and t > s > 0, 

:= Ejf s (Br k! Dr k) , r:(t) ■.= cm/r, 

Cl(s, t) := Cl(t) - Cl(s), Z[( S , t) := T^t) - Z[(s). 

The following property is useful when proving other results of the section. Only minor adjust- 
ments in the proof of \1'6\ Theorem 5.1] are needed to establish it. 

LemmaS. By (A.3) and (A.4), asr^oo, (£ r (-), (xiX(-)), {X2X {■))) => fa(-) *0, p(-),a(-)), 
where Xi( x i,X2) ■= x\, X2{xi, x 2 ) := x 2 , and rj(t) := trj, p(t) := tp, a(t) := to. 

7.2 Compact containment 

The property we prove here, together with the oscillation control result that follows in Sec- 
tion [731 implies tightness of the scaled state descriptors. 

Lemma 9. By (A.3)-(A.5), for any T > and e > 0, there exists a compact set K C M J such 
that 

Hm r .^ 0O P r {if(i) G K for all t G [0,T]} > 1 - e. 

Proof. Fix T and e. It suffices to show that, for each i, there exist a compact set K{ C M such 
that 

lim r _ >00 P r {2[(t) G if* for all t G [0, T]} > 1 — e/I. (29) 
We use the following criterion (see e.g. |15t Theorem 15.7.5]). 

Proposition 1. A set Ai C M is relatively compact if and only if supgg^ < 00 and 

sup ?gA 4 \ [0, n] 2 ) — > as n -> 00. 

Note that 

^(i)(R 2 + ) = %(t) < Z[(0) + ^(T) = ^(0)(R*) +£l(r)(R2_). (30) 

Also note that, if the residual size (patience time) of a flow at time t exceeds n, then its initial 
size (patience time), must have exceeded re, too, which implies the following bound: 

Ximl \ [0,n] 2 ) < 2[(0)(P4 \ [0,n] 2 ) +Z[(T)(R 2 + \ [0,re] 2 ). (31) 
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The sequence {JS^(O) + (T)} re iz converges and hence in tight, i.e. there exists a compact set 
K[ C M such that 

inf r6B P r {I[(0)+£;(T) G K[} > 1 - e/I. (32) 

Put 

K'l := {i G M: for some f G < and 

ml \ [0,n] 2 ) < \ [0,n] 2 ), n G N}. 

Then the criterion of relative compactness for K'l follows from that for K[, and (|30p - (|32p im- 
ply (f29|) with K. L taken as the closure of K'l . □ 



7.3 Asymptotic regularity 

This section contains three Lemmas. Lemmas [TU] and [TT] prove that neither initial nor newly 
arriving flows concentrate along horizontal and vertical lines. These two results are combined in 
Lemma [T2l that implies the oscillation control result of the next section, and also is useful when 
deriving the limiting equations for the state descriptors in Section 17.61 

Recall from Section [5] that, for b > a > 0, 

H b a = R + x[a,b], V* = [a,b]xR+, 
and introduce similar notations 

:= R + x [a, oo), V™ := [a, oo) x R+. 

Lemma 10. By (A. 5) and (A. 6), for any 5 > and e > 0, there exists an a > such that 
hm^nsup^ ||^(0)(^ +a )|| V \\T (0)(^ +a )|| < 6} > 1 - e. 

Proof. Fix 5 and e. Since, for any £ G and a > 0, 

su P , eR+ ||^ +a )|| v \\t(y* +a )\\ < 2su PneN HC^dJII v U(v ( T- l)a )l 

it suffices to find an a such that 

lim^P^TtO) GMJ > 1-e, 
where A4 a := {£ G M 1 : su PneN ||£(# ( n n a _ 1)o )|| V IK(V£i 1)a )|| < 5/2}- 

The set A4 a is open because £ fc A £ G A4 a implies that £ fc G 7W a for k large enough. Indeed, 
pick an N G N such that U(H^ a )\\ V ||£(V^°)|| < 5/2. Then, by the Portmanteau theorem, 

lim^SUp^ ||^(tf(T-l)a)ll V ll^^r-Da)!! 
<E5*-oo m^Kn^ ||e fe (^ ( r-l)JII V ll^(7-l)JII V Wa)ll V ||^W )|| 

<m»K„a lle(^(T-i)a)ll V ||^ ( r-i)a)ll V lk(^ a )|| V IKMfJH < 5/2. 

By (A. 6) and |13] Lemma A.l], there exists an a such that P{£° G M a } > 1 - e. Then, again 
by the Portmanteau theorem, 

lim^^firiO) g Ma} > F{C° eM a }>l-e. □ 

Besides being used in the proof of the fluid limit theorem, the following result is also used 
when establishing convergence of the stationary distributions of the scaled state descriptors, see 
Section [5j 



20 



Lemma 11. By (A. 3) and (A. 4), for any T > 0, 5 > and e > 0, f/iere exists an a > suc/i 
lim^P^sup^ su PxeR+ ||r' new (t)(^ +a )ll V ||^ new (t)(T^+ a )|| < 5} > 1 - e. 

S v ' 

Proof. Fix T, <5 and e. We first construct auxiliary events Jig such that lim r ^ c ^F r > 1 — e, 
and then show that fi£ D for all r, which implies the theorem. 

Definition o/Og. By Lemma there exists a compact set K C M J such that 

lim^oX^^) G ^ for all t G [0,T]} > 1 - e, 

V V ' 

=: Q[ 

and by Proposition [J M := sup^ K U(M. 2 + ) \\ < oo and sup feA - \\£(R%) \ [0,L] 2 || < 5/4 for 
a large enough L. 

For each i, the rate function Aj(-) is positive on {z G z« > 0} and, by Lemma [H it is 
continuous there. Hence, 

A* := mini<i</inf{Ai(^): z { > 5/4, \\z\\ < M} > 0. (33) 

Put 

6 at a 7(A, A 1) 

7:= — —AT and a := . 

72||rj|| 3 

Also pick an N large enough so that 

Na > L + (|H| V 1)T. 

For m, n G N, define the sets 

Im,n '■= \{jn — l)a, ma) x [(n — l)a, na), 

jm,n . = ^ m _ 2 ) + aj ( m + ^ x ^ n _ 2 ) + aj ( n + 1 j j j 

and pick functions <7 mj „ G C K 2__^ 1 j such that 

!/„,„(■) < 9m,n(-) < !/-,«(•)• 
Since is a vector of probability measures, 

£ JK 5m ,n,#}||<|lE pN ^ ,n )ll<9. (34) 

By Lemma[8]and the continuous mapping theorem, for all m, n G N, {g m ,n, v(')(9m,n^) 
as r — >• 00. Since the limits are deterministic, we have convergence in probability. Since the limits 
are continuous, we have uniform convergence on compact sets. Hence, 

lim r ._ >00 P r {maxi< m , in < A rsup ie [o iT ] \\{g m ,n, £ r (t)) -in* {g m ,n,9)\\ < <V(16iV 2 )} = 1. 

V v ' 

=: ni 

Similarly, by (A. 3), 

lim^ooP^sup^^] \\E r (t) - trj\\ < 5/16} = 1. 

* v ' 



For all r, put 

then lim r ^ 00 P r {rig} > 1 — e, and it is left to show that f2g C 



n r := ^nff 2 n^, 
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Proof of fig C f2£. Fix r £ 1Z, t £ [0,T], x £ M+ and i. Also fix an outcome U) £ £1$. All random 
objects in the rest of the proof will be evaluated at this to. We have to check that 

Z[' new (t)(H^ a )<6, (35a) 

Zl' DeW (t)(V* +a ) < 5. (35b) 

We will show (|35ap . ()35bp follows similarly. 

Define the random time r := sup{s < t: Z[' new (s) < 5/4} (supremum over the empty set equals 

y now 

by convention). Although in general r is not a continuity point for (■), we still can 

estimate Z£ (r): 

^I' ncw (r) < 5/2. (36) 

Indeed, if r = 0, then Z^' new (r) = 0, and ([36]) holds. If r > 0, pick a r' £ [(r — 7)+, r] such that 
Z-' ncw (r') < 5/4. Then, by the definition of U r 3 , 

Y t ' ncw (r) < Zl' ncw (r') + (E^t) - < 5/4 + Vi (r - r') + 5/8 < |M| 7 + 35/8, 

and (|36p holds by the choice of 7. 

Now, if t = t, then (|36[) implies (|35ap . and the proof is finished. Assume that r < t. Then, by 
the choice of L and (1361) , 



. new 



(*)(^ +o ) < ?r (t)(m +a n [o,l] 2 ) + 5/4 

E[(r)+1 - 



< 5/2 =: s k 

and in order to have (|35ap . it suffices to show that 

S : = \ T%%>» »» = E m ,„ 6B ; Ej! )+1 **> S «/4. (37) 



— • E mn 

First note that 

S mjn = if m > N orn > N. (38) 

Indeed, consider a flow on route i that arrived at C/^ £ (r, t] with (B^ k ,Dl k ) £ I m . n . If m > TV, 
then > L + ||m||T by the choice of JV, 1?^ — , U 7 ik ,t) > L by the rate constraints, and 
s fe = 0. If n > N, then > L + T by the choice of N, D r ik - (t - U r ik ) > L and again s k = 0. 

Now we estimate E m n for 1 < m, n < N. Fix m, n. Consider two flows k < I such that UT, , £/?J € 
(r,t] and (B[ fe , ^),(^, DJ) £ / m , n . In (r,t], ^(-) > Z[' new (-) > e/4 and ||Z r (-)|| < M, and 
then (|33p implies that 

mi se(T,t] ^(^(s)) > A*. 

If^-C/4> 7 , then 

> 3a < a 

(B£ - S t (Z\ Ul h t)) - {B\ k - S^T, U{ k ,t)) > tA? - (ijT^) > 2a, 
PS - (i - OS) - (14 - (t - U[ k )) > ^ - p[ fc -^) > 2a, 

> 3a < a 

and hence at most one of s k and S; is non-zero. This implies that all arrivals to route i during 
(r, t] that correspond to non-zero summands in S mj „ occur actually during a smaller interval 
(t m ,n,tm,n + 7] f= ( t j*]- Then, by the definition of O^, 

< - 2^ k=Entmn)+1 h m JB lk , D lk ) < S u Ps6[0>T _ 7] - 2^ k=EHs+1) g m , n (B lk , D lk ) 

= su Pse[0,T- 7 ]((fi'm,n,^-(s +7)) - (g m ,n, C\ (s))) < 7^i(5m,n, + 5/(8iV 2 ). 
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We plug the last inequality and (|38|) into £ = Xlm neN ^m,m then (f37|) follows by (f3l|) and the 
choice of 7. □ 

The previous two lemmas are summed up into the following result. 

Lemma 12. By (A.3)-(A.6) ; for any T > 0, 5 > and e > 0, there exists an a > stic/i i/icrf 
lim r ^ oo P r {sup t6[0>T] su PxeR+ \\Z r (t)(H* +a )\\ V ||Z r '(t)(y a f+ a )|| < 5} > 1 - e. 



sup xeR+ \\Z r ^(t){H* x + a )\\ V \\Z r >' mi \t)(V* +a )\\ 



Proof. Note that 

sup :reR+ \\Z r,im \t) 

<su Px6R+ nr(o)(^ +a )n v \\r(o)(v* +a )\\. 

Indeed, 

^' init (t)(i/r a ) < %mH^ +t ) ^d ^ ,init (t)(^ +a ) <^(o)(Ct!#J)' '^- 

Then the lemma follows by ~Z? (■) = (j£ r,vaA _|_ 2 r ' ncw )(.) and Lemmas 1 1 01 and 1 1 1 1 □ 
7.4 Oscillation control 

Here we establish the second key ingredient of tightness of the scaled state descriptors, the first 
one is proven in Section [7,21 

Lemma 13. By (A.3)-(A.6), for any T > 0, 5 > and e > 0, there exists an h > such that 

iim r ^ 00 P r M^ r ', h, T) < 5} > 1 - e, 

V v ' 

where oj(Z t , h, T) := sup{di(Z r (s), Z^ {t)) : s,t € [0, T], \s — t\ < h}. 

Proof Fix T, 5 and e. By (A.3), 

lmWooP^sup^o^] \\E r (t) - tri\\ < 5/4} = 1. 

V ■ v ' 

=: ill 

By Lemma [T21 there exists an a > such that 

fe P _ooP r {sup te[0l ri \\Z r (tW U V a )\\ <5}>l-e. 

V v ' 

=-M r 2 

Pick an /i such that /i(||m|| VI) < 5 A a and < 5/2. We now show that, for all r € TZ, 

D fij n Oj, then the lemma follows. 

Fix r £ 1Z, i and s, t £ [0, T] such that s < t, t — s < h. Also fix an outcome cj E n f^- All 
random objects in the rest of the proof will be evaluated at this co. We have to check that, for 
any non-empty closed Borel subset B C ]R^_, 

Z[(s)(i?)<^)(i? 5 )+5, (39a) 

Z^)(i?)<Zl( S )(i? 5 )+5. (39b) 

First we check (|39ap . Note that it suffices to show 

^(,)( J B)<^[(r)( J B 5 ) + 5, (40) 
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where r := inf{ii G [s,t]: ~Z^(u) = 0} and infimum over the empty set equals t by defini- 
tion. Indeed, if r = t, then (|40p implies (|39ap . If r < t, then by the right-continuity of Zj(-)> 
%(t)(B s ) = zI(t) = 0, and again flU implies (l39aD . 

Now prove (|40p . If r = s, then ()40p holds. Assume that r > s. By the defintion of 

^OO(B) < n [a, oo) 2 ) + 5. (41) 

Since Si(Z r , s, r) < \\m\\h < 5 A a and r — s < h < 5 A a, 

^( S )(5n[a,oo) 2 )<^(T)(i? 5 ), 
which together with (|41|) implies (|40|) . 

It is left to check (I39b1) . Since , s, r) < < 5 and t — s < h < 5, 

-W i {t){B)<Z r i {s){B s ) + {E r i {t)-%{s% 
and (|39bp follows by the definition oi£L\. □ 

7.5 Fluid limits are bounded away from zero 

Rate constraints provide infinite-server-queue lower bounds for bandwidth-sharing networks. 
First we show that properly scaled infinite server queues are bounded away from zero, and then 
the same follows for bandwidth-sharing networks with rate constraints. 

Consider a sequence of infinite server queues marked by r 6 1Z. At t = 0, the queues are empty. 
To the r-th queue, customers arrive according to a counting process A r {-) and have i.i.d. service 
times {-B£}fceN distributed as B r . Let A (•) := A r (-)/r =^> «(•)> where a(t) := ta and a > 0. 
Also let B r => B, where ¥{B > 0} > 0. Denote by Q r (■) the population process of the r-th 
queue and put Q (•) := Q r (-)/r. 

Lemma 14. For any 5 > 0, there exists a C{5) > such that, for any A > 5, 

P r {inf 5 < 4 < A Cf(t) > C{5)} -> 1 as r ^ oo. 

Proof. Let us first explain the result heuristically. Consider the arrivals with long service times, 
i.e. exceeding a b > 0. During (0, 6/2], there are ra¥{B > 6}6/2 such arrivals to the r-th queue. 
They will leave the queue after t = b, and hence, in (b/2,b], the scaled queue length Q (•) 
is bounded from below by aP{B > b}b/2. Similarly, Q (•) > a¥{B > b}b/2 in any interval 
((n- 1)6/2, nb/2], n € N. 

We now proceed more formally. Pick an b € (0, <5) such that 6 is a continuity point for the 
distribution of B, and 

p := F{B > b} > 0. 

Then, asr^oo, 

Pr := P r {B r >b}^p. 
Partition (0, A] into subintervals of length 6/2, 



(°' A ^U 1 < n < mA) ((--W2,n6/2]. 



i<N(A) 

Denote by ~A^ n the scaled number of arrivals during ((n — 1)6/2, n6/2], and by A^ n {b) the scaled 
number of arrivals during ((n — 1)6/2, n6/2] with service times at least 6, 



<:=A r (n6/2)-T((n- 1)6/2), 

-rr . . _ 1 Y ^A-(nfe/2) 

A nW ■--2^ k=Ar{{n _ 1)b/2)+1 hBl>b}- 
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By A =>■ a(-) and p r — > p as r — > oo, 

(AT,...,Z^ (A) ) (ab/2,...,ab/2), 
(A[{b), . . . ,T N(A) Q>)) (apb/2, . . apb/2). 

Pick a C(<5) < apb/2, then 

P r {mf 5 <t< A Q T '(i) > C(S)} 

> F r {mf t€{{n _ 1)b/2M/2] Q r (t) > C(5), n = 2, ... , 7V(A)} 

> P r {X^(fe) > C(<J), n = 1, ... , JV(A) - 1} -»■ 1 as r -> oo. □ 

We can now prove easily that all fluid limits are bounded away from zero outside t = 0. 
Lemma 15. For any 5 > 0, i/iere exists a C{5) > suc/i i/iai, for any fluid limit (Z, Z){-), 

a.s. inf t > 5 mini<j</Zi(i) > C{5). 

Proof. Consider a flow k on route i in the r-th network. By the rate constraints, this flow will stay 
in the network at least for B r ik /mi A D r ik since its arrival. Hence, the route i population process 
Z?(-) is bounded from below by the length Q\{-) of the infinite server queue with arrivals E[ (•) 
and i.i.d. service times {B'- k /m,i A D[ fe }j. e N- Assume that Q[(0) = and put Qj(-) = Q[(-)/ r - 
Then, by Lemma PHI for any 5 > there exists a C(<5) > such that, for any A > 5, 

P r {mf te[5 ,A] mini^^/^Ct) > C(5)} > F r {inf te[(5iA] mini<i</^(t) > C(5)} 1. 

Now consider a fluid limit (Z,Z)(-) along a subsequence {(i? 9 , Z 9 )(-)} gg g. For any compact 
set If C R+, the mapping 99^ : Ds + ^i — > R, ^pk(x) ■= infteK niin g <j</ x(t) is continuous at 
continuous x(-). Hence, f[s^\{Z q ) => f[s.A](Z) and, by the Portmanteau theorem, 

where A > 5 is arbitrary. Then the lemma follows. 

Note also that the constant C(5) does not depend on a particular fluid limit (Z, Z){-). □ 

7.6 Fluid limits as fluid model solutions 

Here we show that fluid limits a.s. satisfy the fluid model equation ([3]). 

Let (Z,Z)(-) be a fluid limit along a subsequence {(Z 9 , Z q )(-)} q( zQ. Lemma [T2l implies that 
(c/. the proof of [111 Lemma 6.2]) 

a.s. Zi(t){d A ) = for all t > 0, alH and A e C, (42) 

where <9^ denotes the boundary of A Then, when proving ([3]) for (Z, Z)(-), it suffices to consider 
sets A from 

C + := {[2, 00) x [y, 00) : x Ay > 0}. 
It also suffices to consider i from a finite interval [0, T]. 

The rest of the proof splits into two parts. First we derive dynamic equations for the prelimiting 
processes (Z q , Z q )(-), and then show that these equations converge to ([3]). 
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Prelimiting equations Fix q € Q, i, t < T and A 6 C + . Fix also an outcome uj € £l q . In 
what follows up to equation ()45|) . all random elements are evaluated at this uj. We have 

Z q (t)(A) =^(0)(^ + (5 4 (Z 9 ,0,t),t)) 

=jSfc 

+ \ EST ^ - ^ ^ -(*-*&))■ (43) 

V v ' 

=: S 

Fix a partition partition < to < *i < ••• < = t, then 

E = 1 yiEtfto) 1 *-^iV-l 

~~ q 2^ fc= i Sfc q 2^/j =0 2^/ k =Ef(t j )+i Sk ' 
Suppose that a function ?/(•) is non-increasing in [to,t] and that, for some 5, 

su Pse[t ,t] \Si(Z q ,s,t) -y(s)\ < 5. 
Now we can estimate S. If U9, £ (tj,tj+i\, then 

Bfk ~ (V<ti) + ^ B ?k ~ S(Z q , Eft, t) < B% - (y(t J+1 ) - 5), 
Dl -(t- t 3 )) < Dl - (t - Eft) < Dl - (i - t j+1 ), 

and 



I) 

'j=o q *-^k=E?(tj)+i 



which can be rewritten as 



„ x-^vAT— 1 „ 



Put 



(44) 



X* := sup AeC su Po < s <,< r \\(C q (s,t)(A) - (t - s) V * 
then, by (l4"iD and (gg]) , 

EjJo 1 " + + ^ - - xq ) 

< Z q {t)(A)-Z q M{A + (S % (Z\^t),t)) (45) 

< m to + X 9 + (^(tj+i - tj)9l(A + (|/(t i+ i) - <5, t - t j+1 ) + x« 
To summarize, we have shown that, for all q € Q and well', 

(Z 9 (-),^)e^, (46) 

where A q C D R+ ^ M / x M + is the set of pairs (£(•), x) such that, for any set A G C + , any 
partition 0<io<^i<---<£iV = £<^ an d any function y(-) that is non-increasing in [to,t] 
and that satisfies sup se r t0)t i ^((l^), s,t) — y(s)\ < 5 for some i and 5, 

E^ {^(tj+i - tj)9 q (A + {y{tj) +S,t- tj)) - x) 
<C(t)(A)-Ci(0)(i4 + (5i«l,C>,0,t),t)) 

< 17^0 + x + E^Tq 1 ~ hWM + - * - <i+i) + a?) • 
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Limiting equations By (A. 3) and (A. 4) (cf. the proof of [ITJ Lemma 5.1]), 

X q =>- as q — >• oo. 



Since the limit of X 9 is deterministic, then the joint convergence (Z q (-),X«) (£(•), 0) holds. 
By the Skorokhod representation theorem, there exist random elements {-2 9 (-)} ge g, an d 

{A" 9 }ggQ defined on a common probability space (fi,.F,P) such that (Z*(-),X*) = (£*(■), -X" 9 ), 

q G Q, and iJ(-) = Z(-), and 

a.s. (Z 9 (-),X 9 ) -»• (Z(-),0) asg^oo. (47) 

Introduce also the total mass processes Z 9 (-) := (1,-Z 9 (-)), g € Q, and Z(-) := (1, £(•)). 
By Lemma [13 (02} and (^6]) . 

a.s. Zj(t) > for all t > and all i, (48a) 
a.s. Ji(t)(5A) = for all t > 0, all i and A £ C, (48b) 
a.s. (i g (-),X 9 ) G .A' for all q £ Q. (48c) 

Denote by Q* the set of outcomes w G O for which (|47p and (|48p hold. We will show that, for 
all uj G f2*, all i, t G [0, T] and ^ G C+, 

ii(t)(A) =^(0)(A + (5i(Z,0,t),t)) 

/■* - (49) 

+ rn / (9 i (>l + (5 i (Z,s,t),t-5))ds, 
J o 

and that will complete the proof of Theorem 

Fix t G [0,T], i and A G C + . Also fix an outcome uj G f2*. All random elements in the rest of the 
proof are evaluated at this uj. 

By (g7D and (l48b]) . 

if (t)(A) -)■ Zi(t)(A) as Q -»■ oo. (50) 
By (|48a[) . the rate constraints and the dominated convergence theorem, 

Si(Z q , s,t) — > Si(Z , s,t) for all s G [0, t] as q — > oo, (51) 

which in particular implies that 

Z?(0)(A + (Si(Z«,0,t),t)) ^Zi(0)(A + (Si(Z,0,t),t)) as ? -> oo. (52) 

Fix ^ G (0, t) and 5 > 0. By ()48a|) . the function 5j(Z, • , t) is continuous in [to, t], and the functions 
Si(Z q ,-,t) are monotone. Then the point-wise convergence (|5ip implies uniform convergence 
in [to,i], and for (/ large enough, 

sup s6 [ t0)t] \Si(Z g ,s,t) - Si(Z,s,t)\ < 5. (53) 

Now fix a partition to < t\ < . . . < = t. The bound (p)3"P and (|48c|) imply that (in the 
definition of A q we take ?/(■) = Si(Z, -,t)) 

< Z?(t)(A) - ZfiA + iS^At)^)) (54) 

< Vi t + x q + J^" 1 (vi(t j+1 - t^A + (S l (Z, t j+1 ,t) -5,t- t j+1 ) + x«) . 

Since 9{{- x R + ) and #i(K+ x •) are probability measures, the set of B G C for which 6i{ds) > is 
at most countable. By (|4*gp . £j(Z, is strictly monotone in [to,*] - Hence, the set P of s G [to>*] 
for which 9i(d A+m2 ,,,t)+S,t-a)) > or ^( a A+(5,(Z, s ,t)-5,t- s )) > is at most countable, too. 
In (|54|). let (/ — > oo assuming that the partition contains no points from T>. Then, by (|47|) . (|50|) 
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and (|52D . 

22 =0 " + *J> *) + 5 ' * " 

< ^(i)(A)-^(0)(^ + (^(Z,0,t),t)) (55) 

< ??i*o + z^ J= o ViW+i ~ tj)9i(A + (Si(Z,tj +1 ,t) -S,t- t j+1 ). 
Now, in ()55p . let the diameter of the partition go to keeping to fixed. Then 

r)i [ 6i{A + (Si(Z, s,t) + S,t- s))ds 
< Zi(t)(A) - Zi(0)(A + (Si(Z, 0, t),t)) 
<Vi t Q + m / Oi(A + (Si(Z,s,t) - 5,t - s)ds. 

Jto 

Finally, in the last inequality, let 8 — > (recall (|48b|) ) and to — * 0, then (|49p follows. 

8 Proof of Theorem [6] 

By the discussion following Theorem [6] and Lemma[2l it is left to show tightness of the scaled sta- 
tionary distributions. It suffices to show coordinate- wise tightness, so fix i. By |14} Theorem 2.1] 
and [El Theorem 15.7.5], the sequence {3^ , jreft is tight if 

sup reTC E r F^ < oo, (56a) 

lim n ^ oo E^(F n oo ) = 0, (56b) 

Wm^WyKH™) = 0, (56c) 

where V£° = [n, oo) x M + and = R + x [n, oo). 

First check (|56ap . For each r, the route % population process Z^(-) is bounded from above by 
the length Q\{-) of the M/G/oo queue with the following parameters: 

(Q.l) at t = 0, there are Z^(0) customers whose service times are patience times of the initial 
flows on route i of the r-th network; 

(Q.2) the input process is the route i input process of the r-th network; 

(Q.3) service times of newly arriving customers are patience times of newly arriving flows 
on route i of the r-th network. 

For all r and t, Z\(t) < Q\{t). As t -> oo, Z\(t) => Y{ and Q\(t) ILffiWD?). Hence, 
Y{ < st n(?7[E r £>[) and E r F- < r/[E r D[/r -> t^EA as r -> oo, which implies (15551) . 

Now check ([56b]) . Note that, if at some point the residual flow size is at least n, then the initial 
flow size was at least n, too. Hence, Z[(-)(V£°) is bounded from above by the length Ql' n (-) 
of the M/G/oo queue whose initial state is as in (Q.l), newly arriving customers are newly 
arriving flows on route i of the r-th network with initial sizes at least n, and service times of 
newly arriving customers are patience times of the corresponding flows. In particular, the input 
process for this queue is Poisson with intensity rjl¥ r {B^ > n}. 

Let fn(') be a continuous function on such that 

Iv^ 1 (-)</n(0<Iv»(-)- 

Then, for all r and t, 

(f n ,Zl(t)) < Z!(t)(V™) <Ql' n (t) 
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Letting t — > oo, we obtain 

XXK+i) < (fn,yD <st m® r {BZ > n}K r Dl), 

and then (|56bp follows. 

Finally, ()56cp is valid due to the following lemma. 
Lemma 16. For any r ElZ, i and Borel set S C R_|_ ; 

yr(R + x S) < st U( V VE r DVW{Dl G 5}), 
where D[ has density P r {£>[ > x}/E r D r i; x>0. 

Proof. Fix r EH, i and a Borel set S C R_j_. It suffices to show that, for any 5 > 0, 

yi(R + x s) < st nfa^rr^r g 5*}), 

so fix 5 > 0. 

Consider the upper bound queue Q\{-) with parameters (Q.1)-(Q.3). Denote by Q r i {t){S 5 ) 
the number of customers in this queue whose residual service times at time t are in S 5 . Then 

Given at time t there are k customers in the queue, denote by D\(t),. . . ,Df.(t) their residual 
service times. By |26^ Chapter 3.2, Theorem 2], 

lim^ 0O P r { J D 1 (f) < Xl ,...,D k (t) <x k \Q r S) = k} = V{Dl < Xl } . . .¥ r {D r t <x k }, 

which together with Ql(t) II(r/[E r D l r ) as t — > oo implies that 

Q:(t)(S s ) ^U^D'^iD: e S s }). 

Let gs be a continuous function on R+ such that 

Ir+xsO) < gs(-) < Ik+xs*(0- 

Then, for any t, 

(g s ,Zl\t)) < Z!(t)(R + x S 8 ) <Q:(t)(S 5 ), 

and as t — y oo , 

^[(R + xS)< ( g & ,yi) < st n«E r Drr{ J Dr G S s }). □ 

Appendix 

Proof of LemmaU^ It suffices to show that, for a vector z G R^_ with the first I' < I coordinates 
positive and the rest of them zero, and a sequence {z fc }fceN C (0, oo) 7 such that z k — > z, we have 
A{z k ) -> A(z). 

Suppose that z fc — > z but A(z fc ) A(z). Since {A(z k )}k^N is a subset of the compact set 
{A G Rj_ : ||A|| < ||C||}, without loss of generality we may assume that A{z k ) — > A' ^ A(z). 

Recall that A(z) is the unique optimal solution to 

maximize ZiUi(A{/ Z{) subject to A A < C, A < m * z, (57) 

where, by convention, Aj/0 := and x (— oo) := 0. 
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For all k, AA(z k ) < C and A(z k ) < m ■ z . Hence, A' is feasible for (|57p and A ■ = = Aj(z) for 
i > I' . Since A' ^ A(z) is not optimal for ([57]) . 

I := Y*' 1 ZiUi(Ai(z)/zi) > V 7 ' i ZiUiik'Jzi) =: r. (58) 

Now we construct a sequence A fc — >■ A(z) such that A k is feasible for the optimization prob- 
lem ([57]) with z fc in place of z. Introduce vectors C k G with = ^i=i'+i ^4^* . 
Put the first /' coordinates of A k to be A k = (Ai(z) - \\C k \\) + A iriizf, and the rest of them 
A^ = Ai(z k ). That is, in the bandwidth allocation A(z), the bandwidth C k , which is required 
for the last I — I' routes, is taken away from the first I' routes. 

Since z k — > z, A k — >■ A(z) and A(z k ) — > A', 

Y!' z k U l (A k /z k )^l and U i {A l {z k )lz\ ) r. 

* — <i=\ — 'i=i 

Also, for all k, 

Y! i=I , +1 4Mtii/zt) = E^, +1 *?z4(a«(**)/*?). 

Then, by for A; big enough, 

i '1 = 1 ^ '4 = 1 

which contradicts to A(z fc ) being optimal for (|57p with z fc in place of z. □ 

Proof of Corollary [7J Fix an FMS (£,z)(-). By Theorem O z(i) — )• z* as i — > oo. Further we 
prove that z(t) — > z* implies £(t) — > (*. It suffices to show that, for any e > 0, there exists a t £ 
such that, for all t > t e , i and Borel sets A C 

Ci(t)(A)<Cf(A e )+e, 

cr(A)<Ci(<)(A 6 )+ e , 1 j 

so fix e > 0. 

For any 5 G (0, mini<j</ z*), there exists a r,5 such that, for all i > t$, 

z * -5:=(z*-6,...,Zj-5)< z(t) < (z* + S, . . . , z? + <5) =: z* + 5. 
Then, for all t > s > r<5 and i, we have 

r;(z* - <5, z* + $)(t - s) < Si(z, s, t) < Ri(z* -5,z* + <5)(t - s), 

^^^^^^^—^ ^>^^^^^^_^ ^^^^^^^^^ 

V V 1 

= : rf =: i?f 

which, when plugged into the shifted fluid model equation (j5a,j) . implies that, for all t > t$, i 
and Borel sets A C 



C*(t)(A) < Ci(ra)(M^ x [i - r tf , oo)) + 77^ / * ^(^ + (rfs, a))ds, (60a) 

J o 

Ci(*)(^) > »7i / 5 ^ + (^fs, s))ds. (60b) 
J o 

Recall from Section that, for all i and Borel sets ^4 C R^, 

/>oo 

= ^ / + (Mz>, (61) 

J 
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Now, there exists a t r e such that, for all i, Borel sets A C ]R^_ and 5 G (0,mini<i</ z*), 

roc roc 

7] { \ 6i(A + (rfs, s))ds < m P{A > s}ds < e/2, (62a) 
Jt' e Jt' e 

poo poo 

th / 6i(A + (Ai(z*)s, a))da < Vi / P{A > s}ds < e/2. (62b) 

Jt' e 

Take 5 € (0, mini<j</ z*) such that 

- \{z*)\\t' e < e/2 and (jr 5 - \{z*)\\t' e < e/2. 
Then, for all i and Borel sets A C P5j_, 

^ / " 0*^4 + (rfs, s))ds <Vi [ e 0i(A £ + (Ai(«*)a, s))ds (63a) 

t?, f " + (Ai(«*)a, s))ds <Vi (' Oi(A £ + (Rfs, s))ds. (63b) 
Also take t" e such that, for all i, 

Ci(Ts)(R 2 + x[^-T 5 ,oo))<e/2. (64) 

Now we put ([6TJ ]) - (j5i|l together in order to obtain ([55]) ; for all t > t £ := (r s + 4) V i", i and Borel 
sets A C R2_, 



j60aT) .(t?>Tt /•t-T« lt62at fti 

Ci(*)(^) < s/2 + rji ei(A + (r 5 lS ,s))d S < e/2 + rji O^A + (rfs, s))ds + e/2 

Jo Jo 

< r/, / 0i(A e + (Ai(z>, + £ < (*(A £ ) + e 
Jo 



and 

l6T1l,ll62at K 



C(^) < »fc / fli(A+(Ai(«>,8))d8 + e/2 
■/ o 

l63bl rt' s l60bl 

< r/i / 0i(A £ + (R!s,s))ds + e/2 < (i(t)(A £ ) + e. □ 

JO 

Proof of Lemma^ For all s < i and e > 0, 

/t pt+x pt+x'+e 

¥{u + x<£<u + x' + e}du = I P{£ > - / IP{^ > u}du 

J S+X J S+x'+e 



CS+X +£ 



< P{£ > u}du < x' - x + e. 



I s+x 



The lemma follows as we first let e — > (applying the dominated convergence theorem) and 
then s — > — oo, t — > oo. □ 
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